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The leader-following consensus of multiple linear time invariant (LTI) systems under switch- 
ing topology is considered. The leader-following consensus problem consists of designing for 
each agent a distributed protocol to make all agents track a leader vehicle, which has the same 
-H 1 LTI dynamics as the agents. The interaction topology describing the information exchange of 
these agents is time-varying. An averaging method is proposed. Unlike the existing results 
■ in the literatures which assume the LTI agents to be neutrally stable, we relax this condition, 
only making assumption that the LTI agents are stablizable and detectable. Observer-based 
leader- following consensus is also considered. 
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1. INTRODUCTION 

Multi-agent system is a hot topic in a variety of research communities, such as robotics, 
sensor networks, artificial intelligence, automatic control and biology. Of particular in- 
terest in this field is the consensus problem, since it lays foundations for many consensus- 
related problem, including formation, flocking and swarming. We refer to survey papers 
[T71 [T§] and references therein for details. 

Integrator and double integrator models are the simplest abstraction, upon which a 
large part of results on consensus of multi-agent systems have been based (see [TH1 [HI 
[U1IE1H1IS]). To deal with more complex models, a number of recent papers are devoted 
to consensus of multiple LTI systems [3ni[23[I3[21[lIJ[Il[Sl[2i[Il[13[lZl[lH]- These 
results keep most of the concepts provided by earlier developments, and provide new 
design and analysis technique, such as LQR approach, low gain approach, approach, 
parametrization and geometric approach, output regulation approach, and homotopy 
based approach. However, most of these results [3UJ ED HI EH E3 M HH H3 mainly 
focus on fixed interaction topology, rather than time- varying topology. How the switches 
of the interaction topology and agent dynamics jointly affect the collective behavior of 
the multi-agent system? Attempts to understand this issue had been hampered by the 
lack of suitable analysis tools. The results of Scardovi et al. [Hj and Ni et al. [TS] 
are mentioned here, because of their contributions to dealing with switching topology in 
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the setup of high-order agent model. However, when dealing with switching topology, 
[2"2"] and [TS] assumed that the system of each agent is neutrally stable; thus it has no 
positive real parts eigenvalues. This assumption was widely assumed in the literatures 
when the interaction topology is fixed or switching. Unfortunately, when the agent 
is stabilizable and detectable rather than neutrally stable, and when the interaction 
topology is switching, there is no result reported in the literature to investigate the 
consensus of these agents. 

To deal with switching graph topology and to remove the neutral stability condition , 
we provide a modified averaging approach, which is motivated by [1] [2j HOj . The averag- 
ing approach was initially proposed by Krylov and Bogoliubov in celestial mechanics |12j , 
and was further developed in the work of [3] [TT] ; for more details refer to the recent book 
[2T1] . Closely related to the averaging theory is the stability of fast time- varying nonau- 
tonomous systems [Tj [TUJ [5] , and more specifically the fast switching systems [33J HI] ■ 
The modified approach in this paper is motivated by the work of Stilwell et al. [23] . 
and also the work of Q] OH [2] ■ Although this work borrows the idea from [23] , the main 
difference of this work from [23 is as follows. The synchronization in [23 is achieved 
under under fast switching condition; that is, synchronization is realized under two time 
scales: a time scale t for the agent dynamics and a time scale for the switching signal 
parameterized by t/e with e small enough. In our paper, we further establish that the 
two time scales can be made to be the same and thus the consensus in our paper is not 
limited to fast switching case. Furthermore, We present an extended averaging approach 
for consensus: a sequence of averaged systems (rather a single averaged systems) are the 
indicators of consensus of the the multi-agent systems. This allows to obtain more relax 
conditions for consensus. At last, We give further investigation on how to render these 
sequence of averaged systems to achieve consensus, and thus ensure the consensus of the 
original multi-agent systems. This was not investigated in [22]. The result in our paper 
shows that if there exists an infinite sequence of uniformly bounded and contiguous time 
intervals such that during each such interval the interaction graph is jointly connected 
and if the dwell time for each subgraph is appropriately small, then consensus can be 
achieved. 

In summary, the contributions of this paper are as follows: 

• Averaging method is applied to leader-following consensus of multiple LTI systems. 

• Results are obtained for a wider class of agent dynamics which is stabilizable and 
detectable than the existing class of neutrally stable agent dynamics. 

• The agent dynamics and the switching signal considered in this paper have the 
same time scale, rather than having different time scales considered in |23j . Thus 
the results in our paper are not limited to fast time switching case. 

The rest of this paper is organized as follows. Section 2 contains the problem formula- 
tion and some preliminary results. Section 3 provides the main result of leader-following 
consensus, and extensions are made in Section 4 which devotes to observer-based proto- 
cols design and analysis. Two illustrated examples are presented in Section 5. Section 
6 is a brief conclusion. 
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2. PROBLEM FORMULATION AND PRELIMINARIES 

This section presents the multi-agent system model, with each agent being a stabilizablc 
LTI system, which includes integrator or double integrator as its special cases. The 
leader-following consensus problem is formulated by use of the graph theory. Some 
supporting lemmas are also included here. 
Consider N agents with the same dynamics 

Xi = Axi + Bui, i = l,2,---,N, (1) 

where Xj € M™ is the agent i's state, and Ui £ M. m is agent i's input through which the 
interactions or coupling between agent i and other agents are realized. The matrix B 
is of full column rank. The state information is transmitted among these agents, and 
the agents together with the transmission channels form a network. We use a directed 
graph Q — (V, £) to describe the topology of this network, where V = {1, 2, • • • , N} is 
the set of nodes representing N agents and £ C V x V is the set of ordered edges 
meaning that agents i can send information to agent j. 
The leader, labeled as i = 0, has linear dynamics as 

x = Ax , (2) 

where xo € K n is the state. Referring agent i g {1, • ■ • , N} as follower agent, the leader's 
dynamics is obviously independent of follower agents. More specifically, the leader just 
sends information to some follower agents, without receiving information from them. 
The interaction structure of the whole agents {0, 1, • • • , N} is described by an extended 
directed graph Q = (V,£), which consists of graph Q, node and directed edges from 
node to its information-sending follower nodes. 

Definition 2.1. (Definitions related to graph) Consider a graph Q = (V, £) with 
V = {0,1,- •• ,N}. 

• The set of neighbors of node i relative to subgraph Q of Q is denoted by Mi = {j € 
V:(j,i)€£,j?i}. 

• A directed path is a sequence of edges {12,13), (^3,24), ■ ■ • in that graph. 

• Node i is reachable to node j if there is a directed path from i to j. 

• The graph Q is called connected if node is reachable to any other node. 
Definition 2.2. (Structure matrices of graph) 

• For a directed graph Q on nodes {1, • • • , N}, it structure is described by its adja- 
cency matrix A G R NxN , whose ij-th entry is 1 if (j, i) is an edge of Q and if it 
is not; or by its Laplacian matrix C = —A + A, where A <G M. NxN is the in-degree 
matrix of Q which is diagonal with z-th diagonal element be | Mi \ , the cardinal of 
Mi, which equals J2j^i a ij- 
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• For a directed graph Q on the node set {0, 1, ■ • • , N}, one uses a matrix H = C + V 
to describe it structure, where C is the Laplacian matrix of its subgraph Q and 
T> = diag^di, ■ ■ ■ , <ijv) with di = 1 if node (0,i) is an edge of graph Q, and with 
dt = otherwise. Obviously, the structure of the graph Q can also be described 
by its Laplacian C. 

It is noted that the graph describing the interaction topology of nodes {0, 1, ■ ■ ■ , N} 
can vary with time. To account this we need to consider all possible graphs {Q p : p G V}, 
where V is an index set for all graphs defined on nodes {0, 1, ■ ■ ■ ,N}. Obviously, V is 
a finite set. We use {Q p : p 6 V} to denote subgraphs defined on vertices {1, • • • , N}. 
The dependence of the graphs upon time can be characterized by a switching law a : 
[0, oo) — > V which is a piece-wise constant and right continuous map; that is, at each 
time t, the underlying graph is Q a (t) ■ 

For each agent i g {1, • • ■ ,N}, if agent j is a neighbor of agent i, the relative in- 
formation feedback to agent i with a gain matrix K to be design later. The 
leader-following consensus problem consists of designing for each agent i € {1, • • • ,N} 
a distributed protocol which is a linear feedback, or a dynamical feedback of 

z t = Y, jeA f iW (xj - Xi) + di{t)(x - xi) (3) 

such that the closed-loop systems CE}-© ac hieve the following collected behaviors: 

]im\\x i {t)-x {t)\\=Q, i=l,---,N. (4) 

t— >oo 

To solve the leader-following consensus problem, the following assumption is proposed 
throughout this paper. 

Assumption 2.3. The pair (^4, B) is stabilizable. 

The following result presents an averaging method for stability of fast time-varying 
linear systems. The general result for nonlinear systems can be found in [TJ [T01 [2] - For 
convenient of our later use, we rewrite the result in the following form. 

Lemma 2.4. Consider a linear time-varying systems x{t) = A(t)x(t) with A(-) : R — > 
R nx ". If there exists an increasing sequence of times tk,k € Z, with tj. — > +oo as 
k — > +oo, tk —> — oo as k — > — oo, and tk+i — tk <T for some T > 0, such that Vt^, the 
following average systems 

x(t) = A k x(t), A k = J -^ ' fc = 0,l,.-. (5) 

i/c+i — ife 

are asymptotically stable, then there exists a* > such that the following fast time- 
varying system 

x(t) = A(at)x(t) (6) 
is asymptotically stable for all a > a* . 
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Remark 2.5. It has been shown in [fj Remark 4] that the value a* can be estimated 
from T by solving the equation 
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for a, where T > is defined above and K v > 0, K > 0, v > are parameters which can 
be determined from the system matrix; furthermore, this equation has for every T > 0, 
K v > 0, K > 0,v > exactly one positive solution a. Now fixing K v > 0,K > 0,v > 0, 
we show that as T — > the corresponding solution a = a(T) — > 0; indeed, T — !> rends 
the right hand side of ([7]) go to infinity, thus requiring — on the left hand side of (J7]) to 
go to infinity, being thus resulting in a — > 0. Therefore, appropriately choosing a small 
T > gives a solution a — a* < 1. 

The following rank property of Kronecker product will be used. The proof is straight- 
forward, being thus omitted. 



Lemma 2.6. For any matrices P,Qi, 
property holds: 
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i/c(P< 



/ Qi 


\ 






( P® 


Qi \ 


Q2 






P® 


Qi 






) 


= rank( 






\ Qn 


J 






\ P® 


Qn J 



The following result will also be used later. 



Lemma 2.7. Consider an n-order differential system x(t) = Aix(t) + A2y(t) with 
Ai e R nxn ,A 2 E W nxm , and y(t) e « m . If A x is Hurwitz and lim^^y^) = 0, 
then limt-i.oo x{t) = 0. 



Proof: Let x(t, xo, y(t)) denote the solution of x(t) — Aix(t) + A2y(t) with initial state 
xn at t ~ 0. Since A\ is Hurwitz, there exist positive number a, 71 and 72 such that 



\\x(t,x ,y{t))\\ < 7ij|a;o||e" 



■l2\\y(t)\ 



where ||y(t)||oo = esssu Pt>o 11^(^)11- Since lim^oo y{t) = 0, then for any s > 0, there 
exists a T > such that 7 2 ||y(£)ll < e/2. Similarly, 7i||xo||e" Q * < e/2. Therefore, 
\\x(t, Xq, y{t))\\ < e. This completes the proof. ■ 



3. LEADERFOLLOWING CONSENSUS OF MULTIPLE LTI SYSTEMS 

This section presents the leader-following consensus of multiple stablizable LTI systems 
under switching topology. Unlike most results in the literature, we do not impose as- 
sumption that A is neutrally stable. For completeness, we first review a result from [TS] 
when the graph is fixed and undirected. 
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Theorem 3.1. For the multi-agent system ([I])-© associated with connected graph Q 
under Assumption 12.31 let P > be a solution to the Riccati inequality 

PA + A T P - 28PBB T P + I n < 0, (8) 

where 5 is the smallest eigenvalue of the structure matrix % of graph (/(which is shown 
to be positive therein), then under the control law Uj = Kzi with K = B T P all the 
agents follow the leader from any initial conditions. 

We now treat the leader-following consensus problem under switching topologies and 
directed graph case. Denoting the state error between the agent i and the leader as 
£j = Xi — xq , then the dynamics of is 

ii = Aei + Bui 

= Asi+BK J2 fa ~ £i) ~ BKdi(t)Ei, i = l,---,N. 
ieJV<(t) 

By introducing e — (e\, e^, ■ ■ ■ , £j[) T , one has 

e = (I N ® A)e - (I N ® B)(C a Q) ® I m )(lN ® K)e - (I N ® B)(V^ t) ® I m )(I N ® K)e 
= [I N ®A-{C a{t) +V a{t) )®{BK)]e 

= [I N ®A-H rj{t) ®(BK)]e. (9) 

The remaining issue is finding conditions on the switching topologies ( i.e., conditions 
on the switching law a) under which one can synthesize a feedback gain matrix K such 
that the zero solution of systems © is asymptotically stable. 

As treated in [5] , consider an infinite sequence of nonempty, bounded and contiguous 
time intervals [tk,tk+i),k = 0,1,-- - , with to = 0, tk+i — £fc < T for some constant 
T > 0. Suppose that in each interval [t/-, tfc+i) there is a sequence of rrik nonoverlapping 
subintervals 

r + l t 2\ r,j \ + m h ,m k + l\ , _ ,1 , _ ,m k + l 

satisfying t]^ 1 — > r, 1 < j < rrik for a given constant r > 0, such that during each of 
such subintervals, the interconnection topology does not change. That is, during each 
time interval [^j^" 1 " 1 ), the graph G a (t) is fixed and we denote it by Qk r The number r 
is usually call the minimal dwell time of the graphs. The r > can be arbitrarily small 
and the existence of such an number ensures that Zero phenomena dose not happen. 
During each time interval [tk,tk+i), some or all of Gkj,j — 1, • • • >mk are permitted to 
be disconnected. We only require the graph to be jointly connected, which is defined as 
follows: 

Definition 3.2. (Joint Connectivity) 

• The union of a collection of graphs is a graph whose vertex and edge sets are the 
unions of the vertex and edge sets of the graphs in the collection. 

• The graphs are said to be jointly connected across the time interval [t, t+T], T > 
if the union of graphs {Q a (s) '■ s £ [t,t + T]} is connected. 
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Assumption 3.3. The graphs Go-It) are jointly connected across each interval [tk, tk+i), k = 
0,1,---, with their length being uniformly up-bounded by a positive number T and 
lower-bounded by a positive number r. 

The following lemma gives a property of jointly connected graphs. When the graph 
is undirected, this result has been reported in [5J H] • We show this result is still valid 
when the graph is directed; its proof is put in the appendix. 

Lemma 3.4. Let matrices W\, ■ ■ ■ , 7i m be associated with the graphs Q\,- ■ ■ , Q m re- 
spectively. If these graphs are jointly connected, then 

(1) all the eigenvalues of Y^Li % have positive real parts. 

(2) all the eigenvalues of Y^iLi T iRi have positive real parts, where > and Y^iLi T i 
1 . 

With this, an averaging method by using Lemma [2^41 is applied to study the stability 
of system whose average system during each time interval [tk, tu+i), k = 0, 1, • • •, is 

x = A k x (10) 

with 

t [In® A — n a(t) ® (BK)]dt 
Afc = — ■ 

tk+l — tk 

= I N ®A-H [tk , tk+l] ®{BK), 

where H [tk , tk+l] = Ete[t»,t fc+ M*) 7 ^)' t j = ~ **)/(*fc+i ~ **)> J =*!.•■ ' >*W- 
Define by Re\ m in{') the least real part of the eigenvalues of a matrix. Define 

8 = mm{ inf ReX min {1lu kM \)\k = 0, 1, • • • }, (11) 

(Tk 1 ,"-iTk mfc )er* 

where where 

ms 

r fe = {(?W , T fe mfc )l X! Tfc 3' T - T J < ^ = ' mfc ^- 
j'=i 

Noting that Re\ rn i n {'fi[t k ,t k+1 )) depends continuously on t^, • ■ • , Tk mk and the set Tk is 
compact, also by referring to Lemma 12.41 one has 

inf Re\ m in{'fi[t k ,t k+1 )) = ReX m in(T k * 1 'Hk 1 H h r^ m ?4 ) > 0, 

[Tfej ,r fcmfc )er fc ■ * - 

which, together with the fact that the set in (jlll) is finite due to finiteness of all graphs, 
implies that 5 in (jTTJ) is a positive number. Then the leader-following consensus control 
can be achieved through the following theorem. 

Theorem 3.5. For the multi-agent system ([TJ-© under Assumption 12.31 associated 
with switched graphs Q a (t) under Assumption 13.31 with T small enough, let P > be a 
solution to the Riccati inequality 

PA + A T P - 25PBB T P + I n < 0, (12) 
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then under the control law Ui = Kzi with K = B T P all the agents follow the leader 
from any initial conditions. 

Proof: We first prove that for each fc = 0, 1, • • • , the average system (fTCTj) is asymptoti- 
cally stable. To this end, let T k e R NxN be an unitary matrix such that T k H [tk . tk+l] T£ 
A , be an upper triangular matrix with the diagonal elements \\ , • • • , Ajy be the eigen- 
values of matrix H[t k ,tu+i\i where T£ denote the Hermitian adjoint of matrix T k . Setting 
x = (Tfe ® I n )x, (US]) becomes 

x = (I N <8> A - A k ® BK)x. (13) 

The stability of (fT3)l is equivalent to stability of its diagonal system 

x = [I N ® A - diag(\\, • • • , \ k N ) <g> BK]x, (14) 

or equivalent to the stability of the following N systems 

Xi — (A — X^BB T P)xi, i=l,---,N. (15) 

Denoting )q = p,* + jv\ , where j 2 = — 1 , then 

P(A - \^BB P ) + {A- \^BB T P)*P 
= P[A {nt + rf)BB T P] + [A- (pi + jvf)BB p ]*P 
= PA + A T P - 2tfPBB T P 

< PA + A T P - 2SPBB T P 

< -I < 

Therefore system (fTU)) is globally asymptotically stable for each k = 0, 1, • • • . 

Using Lemma [2T4l we conclude that there exists a positive a* dependent of T, such 
that Va > a*, the switching system 

i(t) = [I N <8 A - H a{at) g) (BK)]e(t) (16) 

is asymptotically stable. According to Remark 12.51 a* can be made smaller than one if 
we choose T small enough. Since a > a* is arbitrary, just pick a = 1. That is, system 
© is asymptotically stable, which implies that leader-following consensus is achieved. 
■ 

Although the exact value of 5 is hard to obtain, this difficulty can be removed as 
follows. Noting that for two positive parameters 5* < <5, if P > is a solution of 
(p~2|) for parameter 5*, then this P is also a solution of (fT2|) for parameter 5. Thus we 
can compute a positive definite matrix P with a small enough parameter 6* which is 
obviously independent the global information. This treatment has an extra advantage 
that it make consensus control law really distributed since the feedback gain K = B T P 
does not include global information. 

Remark 3.6. During each interval [tk,t k +i), the total dwell time of the m k graphs is 
upper bounded by a positive number T, which is required to be appropriately small to 
make a* < 1. This means that the dwell time of each graph can not exceed a certain 
bound. However, in [15] the dwell time for each graph can be arbitrary since there T is 
not constrained and can be chosen arbitrarily large. 
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Remark 3.7. Note that in (|16l) the switching signal a(at) and state e(t) have different 
time scales, while our result is obtained for system © with a(t) and e(t) have the same 
time scale, and thus the result in our paper is not limited to fast time switching case. 
This distinguishes this work from |23j . 

Remark 3.8. It can be seen that if P > is a solution to (|24l) . then kP,k > 1, is 
also a solution to flU}. Indeed, kPA + kA t _P - 2k 2 6PBB t P + SI n = k(PA + A T P - 
2n5PBB T P + 1/kI„) < k(PA + A T P - 2SPBB T P + I n ) < 0. Therefore, kK is also a 
stabilizing feedback matrix and the n can be understood as the coupling strength. 



4. OBSERVER-BASED LEADER-FOLLOWING CONSENSUS 



This section extends the result in last section to observer-based leader-following con- 
sensus. Consider a multi-agent system consisting of N agents and a leader. The leader 
agent, labeled as i = 0, has linear dynamics as 



x 
Vo 



■ Ax , 
Cx 



(17) 



where yo £ 
as i e {1, • ■ 



is the output of the leader. The dynamics of each follower agent, labeled 
,N},is 



Ui — Cxi 



(18) 



where yi £ M. p is the agent i's observed output information, and Ui £ M. m is agent i's 
input through which the interaction or coupling between other agents is realized. More 
specifically, Ui is a dynamical feedback of Zi. In this section, we assume 

Assumption 4.1. The pair (A,B) is stabilizable, and the pair (A,C) is detectable. 

The observer-based feedback controller is represented as 



Ui 



Ae r 
■■ Fii 



+ K (zi - Zi) + Bu t 



where 



Zt = ^ (Cij ~ Cii) + diCii 



(19) 



(20) 



and the matrices K a and K are to be designed later. 

Remark 4.2. The term in (fT9"]) indicates that the observer receives the output vari- 
able information from this agent's neighbors as input, and the term Zi indicates that 
this observer exchanges its state with its neighboring observers. That is, each observer 
is implemented according to its local sensing resources. Since Zi and Zi are local, the 
observer is essentially distributed, thus then feeding the state of each observer back to 
the corresponding agent is again a distributed control scheme. 
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By further introducing the following stacked vector e = (if , ■ ■ ■ ,eJf) T ,z = (zf, • • • , zf q ) T , 
and by using the structure matrices of graph Q^u), one has 

i = (Jjv ® A)i - [C a{t) ® [K C))i - [D a(t) ® (K C)]e + 

[C a{t) ® (K C)]e + [V a{t) ® (X C)]e + (J* ® B)« 
= [I N (g> (A + BF) - n a{t) ® (X C)]e + ® (K C)][e. (21) 

Then 

s = (I N ® A)e + [I N ® (SJ 1 )^ , 22 x 
i = ® (J^oC)]e + [Jjv ® ^ + Jjv ® (SJ 1 ) - W ff (t) ® (^oC)]e 1 j 

Let e = e — £, that is 




Under this coordinate transformation, system (|22[) becomes 

/ e \ / I N ® (A + 5 J 1 ) I N ®BF We 

U J 1^ Jjv ® A - ?C (t) ® (^C) J ^ e 

Therefore, observer-based leader-following consensus consists in, under jointly con- 
nected graph condition, designing matrices K Q and F such that system ([23]) is asymp- 
totically stable. By separate principle and by referring to Lemma [2. 71 system (1231) can 
by made asymptotically stable through carrying out the following two-procedure design: 

• Design matrix K a such the switched system e = [In ® A — H a (t) ® (-£C)C)]e is 
asymptotically stable; 

• Design matrix J 1 such the i = [In ® (A + BF)]e is asymptotically stable; 

The first step can be realized by referring Theorem 13.51 replacing the pair (A,B) with 
(A T , C T ). The second step is a standard state feedback control problem. We summarize 
above analysis in the following theorem. The rest of its proof is essentially similar to 
that of Theorem [331 an d is thus omitted for save of space. 

Theorem 4.3. Consider the multi-agent systems (|17H18p associated with switching 
graphs Ga(t) under the Assumptions 13.31 14.11 with T small enough. Let P > be a 
solution to the Riccati inequality 

PA T + AP- 2SPC T CP + I n < 0, (24) 




then under the control law (flUf with K Q = PG T and F being such that A + BF is 
Hurwitz, all the agents follow the leader from any initial conditions. 
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5. SIMULATION RESULTS 

In this section, we give two examples to illustrate the validity of the results. Consider a 
multi-agent system consisting of a leader and four agents. Assume the system matrices 



0.5548 -0.5397 -0.0757 \ / 3 5 

0.3279 -0.0678 -0.4495 1,5=1 3 -2 | ,G 
-0.0956 -0.6640 0.0130 



1 -1 2 

-4 2 -3 



We suppose that possible interaction graphs are {G±, G2, G3, G4, G5, Gq\ which are 
shown in Figure [TJ and the interaction graphs are switching as G\ — »■ G2 — > G3 — > 
G4 — » G5 — » Ge — i> Gi —>■••• , and each graph is active for 1/2 second. Since the graphs 
Gi U G2 U G3 and G4 U G5 U Ge are connected, we can choose tk = k, tk+i = k + 3/2 
and t° k = k,t\ = k + 1/2, t\ = k + 2/2, f| = + 3/2 with k = 0, 1 • • • . We choose a small 
parameter (5 = l/3min(0.3820, 0.1732) = 0.0057. The matrices K in Theorem I3~5l and 
Ko , F Theorem 14.31 are calculated as 

0.7520 5.9852 -2.7041 
12.6966 -3.8441 1.6419 



K = 

and 

r , 0.6338 -0.5087 0.3731 . . )r ., 

F={ -0.9077 0.4509 -0.1938 1 < *° = 1 " 9 ' 6111 4 ' 1353 



6.7092 -9.1532 
9.6111 4.1353 
7.7514 -1.1756 



respectively. With the same initial condition, the simulation results of Theorem 2 and 
Theorem 4 are shown in Figure and Figure O respectively. 



6. CONCLUSIONS 

This paper presents an averaging approach to leader-following consensus problem of 
multi-agent with linear dynamics, without imposing neutrally stable condition on agent 
dynamics. The interaction topology is switching. The proposed protocols force the 
follower agents to follow the independent leader trajectory. The result is extended to 
observer-based protocols design. Such design can be separated as a two-step procedure: 
design asymptotically stable distributed observers and asymptotically stable observer- 
state-feedback protocols. 



APPENDIX 

The appendix is devoted to the proof of Lemma 13.41 To this end, we first cite the 
following result. 

Lemma 6.1. (Gersgorin)[7| For any matrix G = [gij\ € M. NxN , all the eigenvalues of G 
are located in the union of N Gersgorin discs 

Ger(G) := U^{z G G : \z - g u \ < £ \ 9ij \}. 
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The definition of weighted graph will also be used in the proof of what follows. If 
we assign each edge of graph Q a weight toy, we obtain a weighted graph Qyj = 
(V, £, W), where W = [it%-]. For an graph £/ and any positive number k > 0, the graph 
/cC? is defined to be a weighted graph obtained from Q by assigning a weight A: to each 
existing edge of Q. For two graphs Q\ and ^2, their union is a weighted graph and the 
weight for edge is the sum of weights for the two edges in the graphs Q\ and 



Q2 respectively. The weighted Laplacian of graph Q y 



is defined as C w = 



-A, 



A, 



where A v 



Q w is defined as % = C. 



% 3 l J 



and A w (i, i) = w ij a iji the weighted structure matrix of graph 



..• .... — «...,. <- w , where C w is the weighted Laplacian of the subgraph 

Q w of Q w , and V w = diag (w 01 di, ■ ■ ■ ,w 0N d N ). 

Proof of Lemma 13. 4t (1) Denote the Laplacian matrix and the structure matrix 
of graph Qi by Ci and T-Li respectively, and denote the Laplacian matrix of graph Qi by 
C. 

We first prove the case when m = 1. By definitions, it can be easily verified the 
following relationship 






••■ 


-di 






Hi 


-d N 





Since the graph Qi is connected, then rank(Ci) = N [18] 
formed by the last N rows of C\ has rank N. Note that 




= T>i 



I 1 



V 1 



£1 



/ 1 



V 1 



■T>i 



( 1 



V 1 



Thus the sub- matrix M.\ 
I 1 



V 1 



that is, the first column of matrix M.\ is a linear combination of its last N columns. 
Therefore, rank{%i) = N . Furthermore, we claim the eigenvalues of the matrix T-Li are 
located in closed-right half plan; indeed, from Gersgorin theorem, all the eigenvalues of 
H are located in 

Ger(H) = u£ x {z e C : \z - l u - d,-| < |JV t -|} , 

and therefore they are located in the closed-right half plan by noting that la = \Ni\ and 
di > 0. We thus conclude that all the eigenvalues of Hi have positive real parts. 

We proceed to prove the case when m > 1. Obviously, for a union Q v of a group 
of weighted graphs {Q%,--- , ^ m }, its weighted Laplacian matrix C v is the sum of the 
Laplacian matrices {£1, • ■ ■ , C m } of graphs {Q\, • • ■ , G m }, and 



/ 






• 













Hi+ ■ 


• +H m 


-< 
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where 



/ 


-d? \ 




( ~d\ \ 




( -d? \ 


V 




- 


\ J 


+ + 


V - d N J 



with (d{,d,2, ••• , d 3 N ) T be the diagonal elements of matrix Vj, j = 1, • • • , m. When the 
graphs are jointly connected, that is, when the Q v is connected, the matrix C u has a 
simple zero eigenvalue. Argue in a manner similar to that of m = 1 case, it can be 
shown that the all the eigenvalues of the matrix ri\ + ■ ■ ■ + V. m have positive real parts. 

(2) Similar discussion as given in (1) for the weighted graphs t\Q\, ■ ■ ■ ,tnQn yields 
the conclusion. 
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Fig. 1. Six possible interaction topologies between the leader and the 

agents. 
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Fig. 2. Simulation for Theorem I3.5I The error trajectories between 
the leader and each agent. 
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(a) Agent 1 follows the leader 
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(c) Agent 3 follows the leader 
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(d) Agent 4 follows the leader 
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Fig. 3. Simulation for Theorem I4.3I The error trajectories between 
the leader and each agent. 



